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Abstract. We apply the Matrix Element Method (MEM) to mass determination of squark pair 
production with direct decay to quarks and LSP at the LHC, showing that simultaneous mass 
determination of squarks and LSP is possible. We furthermore propose methods for inclusion of 
QCD radiation effects in the MEM. 
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The goal of the LHC at CERN, scheduled to start this year, is to discover new physics 
through deviations from the Standard Model (SM) predictions. After discovery of devi- 
ations from the SM, the next step will be classification of the new physics. An important 
first goal in this process will be establishing a mass spectrum of the new particles. One of 
the most challenging scenarios is pair-production of new particles which decay to invisi- 
ble massive particles, giving missing energy signals. Many methods have been proposed 
for mass determination in such scenarios (for a recent list of references, see e.g. [1]). 

In this proceeding, we report the first steps in applying the Matrix Element Method 
(MEM) in the context of supersymmetric scenarios giving missing energy signals. After 
a quick review of the MEM, we will focus on squark pair production, a process where 
other mass determination techniques have difficulties to simultaneously determine the 
LSP and squark masses. Finally, we will introduce methods to extend the range of 
validity of the MEM, by taking into account initial state radiation (ISR) in the method. 

The Matrix Element Method is a procedure to measure a set of theoretical parameters 
from a sample of experimental events. It associates to each event x a weight P(x\oc), 
which is a measure of the probability for an event x to be observed given a set of 
parameters a. The computation of these weights is done by convoluting the production 
of a parton-level configuration y, given by the squared matrix element \M a \ 2 (y), with 
the probability that the parton-level configuration y evolves into the reconstructed event 
x 9 as modeled by a transfer function W(x,y). As a result, the weight of a specific event x 
can be written as 

P {x\a) = — f #(j)|M a | 2 (j)^ 1 ^ 2 / 1 ^ 1 )/ 2 (^)W(x,j) (1) 

where f\(q\) and f^qi) represent the parton distribution functions, d(j)(y) the phase 
space measure and a a the total cross section. The normalisation by the cross-section 
ensures that P(x\a) is a probability density: / P(x\a)dx = 1. 
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In principle, this definition provides the best possible discriminator on a event-by- 
event basis, and can be used to measure a. With a sample of N events, we can construct 
a log-likelihood: 



where the acceptance term Acc(x), corrects for the bias introduced by detector accep- 
tance and event selection.The set of parameters a maximizing the likelihood correspond 
to the most probable value. 

Although the MEM is conceptually very simple, the numerical evaluation of the 
weights is not straightforward due to the large variations of the integrand. Indeed, both 
the square matrix-element \M a \ 2 and the transfer function W(x,y) are highly non uni- 
form on the phase-space. In consequence, a specific Monte Carlo phase-space integrator, 
designed to deal with the specific behavior of the integrand is required. 

MadWeight [2] is a publicly available phase space generator dedicated to perform this 
type of integration. In order to find the best phase-space mapping, MadWeight applies 
a series of local changes of variables to promote, when possible, the invariant mass of 
resonant propagators as variables of integrations. In this way, MadWeight is able to find 
the best phase-space mapping for any decay chain and a large class of transfer functions. 



The MEM has been successfully applied for determination of the top quark mass at the 
Tevatron [3]. However, applying the MEM to supersymmetric scenarios brings new dif- 
ficulties, including a large parameter space, a priori unknown order of particles in decay 
chains, and unknown masses of invisible particles. In the preliminary study presented 
here, we focus on the latter of these difficulties by studying squark pair production, 
where the squarks decay directly to quark and LSP, the lightest neutralino. The signature 
is two hard jets and large missing transverse energy. This scenario is interesting, because 
the mt2 "kink" method [4], as well as other methods for simultaneous reconstruction of 
the squark and neutralino masses often fails for two-body decays of the squarks. 

As a first check that the Matrix Element method can indeed reconstruct masses in 
this scenario, we simulate squark pair production at the LHC, and study this process at 
parton level without cuts, with rriq — 561 GeV and m^o = 97 GeV. The resulting negative 
logarithmic likelihood function is shown for 100 events in Fig. la). The minimum of 
— ln(L) forms a "valley" in the (m^m^o) plane. This valley is closely aligned with 
the expression for the energy of the quark in the rest frame of the decaying squark, 
i.e., the maximum pr of the quark produced by squark production at rest, given by 
p™ x = (m~ — m 2 ~ )/2mq. Fig. 1 b) and c) show the extent of the valley using 3000 and 

7500 events, respectively (corresponding to an integrated luminocity of 10 fb _1 and 25 
fb _1 respectively, assuming 100% branching ratio). With 3000 events, only an upper 
limit on m^o can be determined. For 7500 events, however, the mass range for the %° 
is already quite restricted, and is determined to be between 65 GeV and 160 GeV. The 
/?5? ax value is very well determined, with an uncertainty of ±2 GeV (less than 1%). 
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FIGURE 1. Matrix element - log(L) for a) 100, b) 3000 and c) 7500 parton level qq -> qqX°X° events 
at the LHC. The true values are at (m^, ra^o ) = (561 , 97) GeV. b) and c) show the valley in the (m^o , pf* x ) 
plane. The blue crosses show parameter points with — log(L) at most 4 above — log(L) m i n (corresponding 
to approximately 90% CL). 



Further studies indicate that these results are still valid when realistic cuts, and the 
effect of hadronization and detector simulation, are considered, although the efficiency 
is then reduced. See [5]. 

INITIAL STATE RADIATION 

One challenge associated with the Matrix Element Method is how to deal with initial 
state radiation (ISR). Energetic jets radiated from the intial state partons will be much 
more ubiquitous for high-/^ physics at the LHC than at the Tevatron [6]. To simply 
ignore their effect leads to bad or impossible fits of the matrix element algorithm; on the 
other hand, a very strict veto on additional jets (besides jets expected from the primary 
hard process) would cause a significant reduction of signal statistics. 

Conceptually, the most straightforward way of including events with sizeable ISR 
would be the use of matrix elements with additional partons in the final state. In practice, 
however, there is a limit to the multiplicity of the matrix elements due to the computing 
time for the matrix element generation and the phase- space integration for the weight 
computation. Therefore, in most cases, only one or two jets from radiation can be 
included. 

This situation is illustrated in Fig. 2 for the example of 1000 simulated events of 
di-lepton top pair production at the LHC, pp —> tt + nj —> bbl + l~ViVi + n j, where n j 
indicates n jets from ISR. The events have been generated with Pythia 6.4/PGS4 [7, 8]. 
Backgrounds have not been included. As is evident from the figure, the inclusion of 
matrix elements with one extra jet already significantly improves the agreement (peak 
of likelihood curve) with the true input value m t = 175 GeV and also the statistical 
precision (width of likelihood curve). 

We also propose a different method to treat ISR. We observe that the most significant 
effect of ISR is on the kinematics of the events, since without proper inclusion of ISR 
the momentum balance would be violated. This effect can be taken into account by 
simply boosting the hard event by the momenta of the initial state jets.The longitudinal 
incoming momenta are integrated over in the computation of the likelihood (see (1)), 
so it is sufficient to perform the boost for the transverse coordinates only. Since the 
identification of the b jets from top decay is not unique, we sum over all permutations of 
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FIGURE 2. Reconstruction of the top quark mass from a matrix element likelihood fit to 1000 di-lepton 
top-antitop events at the LHC with y/s = 14 TeV. A top mass of m t = \15 GeV has been used for the event 
generation. Blue: "ideal" case without ISR in the event generation. .Red solid (rightmost curve): realistic 
event generation with ISR, but events with extra jets with pr > 40 GeV have been vetoed. Red dashed: 
same as above but also including tt matrix elements with one additional jet. Green: boost method to deal 
with ISR in the events. 



the four hardest jets in the event. 

The result of a fit with this boost method is shown by the green curve in Fig. 2. It 
agrees with the true input value within statistical errors. The quality of the fit result is 
comparable to the computation with explicit tt + j matrix element (red-dashed curve), 
but requires less computation time. Furthermore, the result of the boost method can be 
improved by including QCD splitting functions Sudakov form factors for the (cumula- 
tive) ISR on each leg [5]. 

In summary, we have here presented a first study in applying the Matrix Element 
Method to supersymmetric particle production with decay to massive invisible particles. 
In particular, we showed that it is possible to simultaneously determine the squark 
and LSP masses in squark pair production with exclusive decay to quark and LSP. 
We also proposed two methods to include the effects of initial state radiation in MEM 
calculations. For more details we refer to [5]. 
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